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ABSTRACT 

Motivated by kinematic observations of young embedded clusters, this paper explores possible kine- 
matic signatures produced by asphericity and departures from initial virial equilibrium in these sys- 
tems. Specifically, the kinematic quantity that is measured and calculated in this study is the distri- 
bution of the line-of-sight velocities as a function of position along the cluster. Although clusters are 
found within a wide range of sizes, we focus on the regime with stellar membership N ~ 10'^. The 
gravitational potential of these young clusters is dominated by the gas, and the geometry of the gas 
distribution is generalized to include axisymmetric (and triaxial) forms. With this loss of symmetry, 
the kinematic results thus depend on viewing angle. This work also considers the stars to begin their 
trajectories with subvirial speeds, as indicated by observations of core motions in such clusters. Our 
results determine the conditions necessary for the kinematic signature to display interesting structure, 
i.e., a non-spherical potential, a viewing angle that is not along one of the principal axes, and subvirial 
starting conditions. We characterize the effects on this signature due to projection angle, initial stellar 
velocities, cluster elongation, and star formation efficiency. Finally, we compare our theoretical results 
to recent kinematic observations of the Orion Nebula Cluster; we find that the observations can be 
explained provided that the cluster is non-spherical, starts with subvirial initial velocities, and is not 
viewed along a principal axis. 

Subject headings: open clusters and associations: general, stellar dynamics, methods: n-body simula- 
tions, stars: formation 



1. INTRODUCTION 

Most sta rs are thought to form within clusters of some 
type (e.g.. iLada fc Ladal [200l iPorras et"an l2003f ). aL 
though a great deal of controversy remains concerning 
the distribution of cluster properties. Nonetheless, given 
that stars form in clusters, two overarching questions im- 
mediately arise. The first considers the clusters as astro- 
physical objects: [1] How can molecular clouds produce 
aggregates of iV > 100 stars with centrally concentrated 
surface density, with the massive stars near the center, 
and with a stellar mass distribution that follows the IMF, 
all within a few pc and within a few Myr? A second vi- 
tal question then becomes: [2] If stars form in clusters, 
how does the cluster environment affect star formation 
and the accompanying process of planet formation? A 
complete understanding of star and planet formation re- 
quires detailed answers to both questions. However, this 
paper will focus on one specific issue within this larger 
context. Observed young embedded stellar clusters dis- 
play departures from both spherical symmetry and initial 
virial equilibrium, and this paper will explore the effects 
of these complications on the kinematics of cluster mem- 



bers through iV-body simulations. These results, in turn, 
will help provide a contribution to the overarching ques- 
tions posed above. 

Departures from spherical symmetry in star forming 
regions have long been known, but have been little stud- 
ied. On the scale of cluster themselves (with radii of 
a few pc, and masses of 100 - 3000 Mp)), for example, 
some systems are categ orized as irregular (jLada fc Ladal 
l2003f) . or filamentary (IWalsh et al.ll2004D. or otherwise 
described as non-spherical ([Alien et al.ll2007l ). Further- 
more, the star formation efficiencies are low, so that the 
gravitational potential in these systems is dominated by 
the mass of the gas (at least in the early phases of evolu- 
tion, before gas expulsion). As a result, the gas potential 
must be generalized to inclu de the observed de partures 
from spherical symmetry fsee lAdams et al.ll2007l and §2). 
For completeness, we also note that on the smaller scale 
of cores (or clumps or kernels) that represent individual 
star formation events, the pre-coUapse gas is observed to 
be non-spherical, with typical aspect ratios of 2:1 (e.g., 
iMvers et al.l[l99ll : iRvdenlflQOl . 

Observations are starting to show that clusters be- 
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gin their evolution with subvirial initial conditions, and 
hence this complication must be included in our simu- 
lated clusters. In many regions, pre-stellar clumps are 
obser ved to move subsonically before the clumps form 
stars jPcrctto ct al."2006t IWalsh et al.l[2003 : lAndrel[200l : 
iKirk ct al, 2007,,). implying that newly formed stars be- 
gin their dynamical evolution with subvirial speeds. Mo- 
tivated by such observations, our work considers clusters 
that are seeded with subvirial stellar velocities. As such, 
our N-body simulations (sec §2) differ f rom those of many 
preceding studies (thoug h not all, see lBonnell fc DaviesI 
119981 lAdams et alJ 12006( 1 that assume the initial phase 
space variables of the stars are close to virial equilibrium. 
Subvirial initial conditions can have a significant im pact 
on the early cluster evolution () Adams et al.l [20061 ). and 
arc thus considered herein. 

The theoretical work presented here provides a deter- 
mination of the kinematic signature of young embedded 
clusters with both subvirial starting conditions and non- 
spherical gas potentials. Fortunately, astronomical ob- 
servations arc now becoming available to compare with 
these theoretical calculations. In this paper, we focus 
on clusters of moderately large size, with stellar mem- 
bership TV ^ lO'^, appropriate for the Orion Nebula 
Cluster (ONC). Clusters of this size are large enough 
for interesting kinematic signatures to arise in observa- 
tions, and small enough that many such clusters are ex- 
pected within the galactic cluster population. In addi- 
tion, the ONC is an example of a nonsphcrical cluster 
as the stellar population is elongated North to South 
along the filament of dense mo lecular gas in the region 
(|Hillenbrand fc Hartmannlll99 8') . Recent observations in 
the ONC (|Furesz et al.l 120081 Tobin et al. 2009, sub- 
mitted to ApJ, T09 hereafter) display a gradient in the 
radial velocity structure along the length of the cluster, 
much like that expected for the elongated subvirial clus- 
ters considered herein. 

This paper is organized as follows. In §2, we outline 
the theoretical approach used in this paper. Specifically, 
we describe the A^-body codes, the required number of 
realizations of the numerical experiments, the implemen- 
tation of subvirial starting speeds, the inclusion axisym- 
metric and triaxial gas distributions, as well as the spec- 
ification of simulation parameters. We then present the 
theoretical results of our simulations in §3, with a fo- 
cus on the velocity signatures produced by the depar- 
tures from equilibrium and spherical symmetry. In §4, 
we compare our results to observations, primarily the 
kinematic velocity measurements recently carried out in 
the ONC. We find good qualitative agreement, and rea- 
sonable quantitative agreement, and suggest that the 
observed kinematic signature requires subvirial starting 
conditions, non-spherical potentials, and viewing angles 
that do not coincide with the principal axes of the sys- 
tems. These conclusions and other results are summa- 
rized in §5, along with a discussion of their implications. 

2. METHODS 

A^-body simulation techniques are used to calculate the 
dynamics of young stellar clusters during the embedded 
star forming epoch. We consider clusters with nonsphcr- 
ical geometries and subvirial initial velocities and study 
the observed kinematics of these clusters. In this section, 
we discuss the A^-body code used, the simulation pa- 



rameters chosen, and the experiments performed to iden- 
tify and characterize kinematic signature in nonsphcrical 
clusters with subvirial initial conditions. 

2.1. N-Body Simulation Techniques 

The dynamical evolution of a young embedded cluster 
depends on its initial stellar distribution, the distribu- 
tion and removal mechanism of the embedding gas, the 
star formation history, and is es pecially sensitive to the 
initial velocity distribution (e.g.. lAdams et a"ni2006f ). In 
this paper we complete a suite of A^-body simulations of 
young embedded clusters to understand how the cluster's 
initial spatial and velocity distribution imprints itself on 
the evolved cluster's kinematic struct ure. The NB0D Y2 
direct integration code developed by lAarsethl (j2001[ ) is 
used as a starting point to calculate the cluster's dynam- 
ics from the star-forming epoch out to ages of 10 Myr. 
As outlined below, we modify the code to include specific 
stellar and gas distributions and star formation epochs 
that are similar to those observed in young embedded 
clusters. 

Distributions of stellar positions and velocities are 
required to characterize kinematic signatures in these 
young clusters. We find that ~ 50 equivalent realiza- 
tions of each set of cluster initial conditions are required 
to provide statistically robust distributions of stellar po- 
sitions and radial velocities. This study focuses on the 
kinematic signatures produced by two different contribu- 
tions to the initial conditions commonly present in young 
embedded clusters: asphcricity (combined with effects 
due to projection of a 3D cluster onto the 2D sky), and 
subvirial initial velocities. In addition, we consider the 
effect of extinction on the observed kinematic structure 
of young clusters. 

Cluster evolution is numerically integrated from the 
star formation epoch (0 — 1 Myr) through the gas re- 
moval phase (at i = 5 Myr) out to ages of 10 Myr (see 
^2.2l for further detail). After 50 realizations of each set of 
initial conditions are completed, the results are combined 
to produce distributions of stellar positions and veloci- 
ties. The simulations provide three dimensional position 
and velocity information for each stellar member at in- 
tervals of 0.1 Myr. These six phase space coordinates 
are reduced to two position coordinates in the plane of 
the sky and two velocity components, one along the line 
of sight (radial velocity) and one in the plane of the sky 
(transverse velocity). For the sake of definiteness, the z' 
projected axis is defined by the projection of the three 
dimensional z axis (the major axis in elongated clusters) 
onto the plane perpendicular to the line of sight. The 
terms "North" and "South" with reference to projections 
of simulated clusters correspond to positive and negative 
z' values respectively. This nomenclature is arbitrary 
and has been chosen to coincide with that of the ONC, 
which is elongated North to South and displays North- 
South asymmetry in the radial velocity structure (see 

m- 

2.2. Simulation Parameters 

Cluster Membership N : We consider moderately large 
young clusters comparable in size to the Orion star form- 
ing region. Specifically, our simulated clusters contain 
A^ = 2000 or 2700 stars. Estimates of the stellar popu- 
lation of the ONC vary depending on the cluster radius 
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adopted, and our choic es of N roughly reflect th e range 
of this variation (e.g. iHillenbrand fc HartmannI (|1998l ) 
advocate N = 2300, near the center of this range). We 
found that there were no significant differences in the 
kinematic signature observed in clusters of 2000 and 2700 
stars. Another motivation for our choice of cluster size 
is that clusters of this size have large enough member- 
ships that kinematic signatures may be identified in the 
data from observed clusters and compared to the results 
presented in this study. 

The stellar masses ar e sampled from the log- normal 
analytic fit present ed by I Adams fc Fatuzzol (|1996( l to the 
standard IMF of iMiller fc Scald ()1979f ). The masses 
range from 0.07 to IOMq. A limited number of simu- 
lations with equal stellar masses were completed and we 
found that the kinematic signature discussed herein was 
present in both the single stellar mass clusters and those 
with a more realistic stellar mass distribution. 

Stellar Distribution: One goal of this study is to com- 
pare kinematical signatures of clusters with spherical ge- 
ometries to those of elongated clusters. Spherical clus- 
ters have centrally concentrated stellar distributions de- 
scribed by a r~^ density profile. Two elongated stel- 
lar distributions arc considered: uniform density pro- 
late spheroids and centrally concentrated clusters with 
p ^ m~^, where = jo} -|- -I- jt? is the 

triaxial coordinate. We compare elongated clusters with 
different aspect ratios to characterize the effect of cluster 
elongation on the observed kinematic signature. Table [1] 
summarizes the initial stellar density distributions used 
in the cluster simulations. 

Initial Speeds: Initial stellar velocities are sampled 
from a uniform distribution within a unit sphere, produc- 
ing an isotropic and position-independent velocity dis- 
tribution. The velocities are then scaled to produce a 
cluster in a particular virial state, defined by the virial 
ratio Q = \K/W\ where K is the total kinetic energy 
of the cluster and W is the total virial potential energy 
of the cluster. A cluster in virial equilibrium is defined 
to have Q = 0.5. Simulations of clusters with virializcd 
initial velocities are compared to clusters with subvirial 
initial velocities. The subvirial clusters are chosen to 
have Q = 0.04 or 0.15, which corresponds to initial ve- 
locities that are approximately one third or one half of 
the virial velocity, respectively. 

Star Formation History: The simulated clusters have a 
star formation epoch that lasts for the first 1 Myr of the 
cluster's evolution. During simulation initialization each 
star, regardless of mass or initial position within the clus- 
ter, is assigned a random formation time between and 1 
Myr (with a uniform distribution, independent of stellar 
mass, over this interval). A star is then tied to its for- 
mation site (chosen as described above) until its collapse 
phase of star formation is complete. That is, the stars are 
initially included in the simulation as point masses, but 
are held at a fixed position until their formation time. Af- 
ter that time, the star becomes free to move through the 
total gravitational potential of the cluster with an initial 
velocity sampled from the distribution described above. 
As a result, the stars do not execute ballistic orbits until 
after they have formed. For completeness, we note that 
stellar evolution is not included in these simulations. 

We vary the star formation efficiency (SFE) from 17% 
to 50%, which spans the range of mass estimates of the 



gas in the region of the ONC. It is important to note 
that although the evolution of a cluster subsequent to 
gas removal depends strongly on the SFE, its evolution 
prior to gas removal is more sensitive to the virial ratio 
than it is to the SFE alone. In other words, provided that 
the stars are moving with sufficient velocities to account 
for the additional potential due to excess gas (i.e., their 
virial ratios are comparable), the evolution of a cluster 
with an SFE of 33% or 50% will be qualitatively similar 
during the embedded stage. Increasing the SFE will, 
however, decrease the average virial velocity. After gas 
is removed from the cluster, the subsequent dynamical 
evolution does depend on the SFE. Of course, increasing 
the SFE will increase the rate of close encounters, but 
this e ffect is modest on t he (short) time scales of interest 
here (| Adams et al.ll2006| ). 

Embedding Gas Distributions: We assume the distribu- 
tions of stars in a cluster roughly traces the geometry and 
density of the embedding gas. T hus a spherical c entrally 
concentrated Hcrnquist profile (jHernquistI 1 1 990( 1 is cho- 
sen to represent the em bedding gas in a spherical cluster 
see lAdams et ai]|2006[) . 

Likewise, the elongated stellar distributions are embed- 
ded in elongated gas potentials. Specifically, the uniform 
density stellar distributions are embedded in a uniform 
density gas distribution that is twice the extent of the 
stellar distribution. The homoeoid theorem states that 
the net force on a particle within a uniform density ho- 
moeoid shell is zero. Thus the larger gas distribution 
allows for a simpler computation of the force and poten- 
tial terms due to the embe dding gas without changin g 
the dynamics of the system (iBinnev fc Tremaindll9 87'). 

The centrally concentrated prolate clusters are embed- 
ded in a static gas potential of the form p ~ m~^, where 
m is the generalized coordinate. Calculation of the force 
terms, analytic expressions for the potential, and a dis- 
cussion of orbits and orbit i nstabilities wi t hin th is triaxial 
potential were presented in lAdams et al.l (|2007f ). The gas 
distributions and associated parameters are summarized 
in Tabled 

Observations indicate that embedding gas does not re- 
main in young clusters for a long time. After a few mil- 
lion years, winds from hot young stars begin to carve out 
the embedding gas and very few embedded clusters are 
found with ages greater than ~ 5 Myr. Our simulations 
account for gas removal as a temporal step function in 
the evolution of the static gas potential: at t = 5 Myr, 
the gas potential is completely removed from the cluster 
which then continues to evolve due to interactions be- 
tween the stars. After gas removal, the cluster expands 
and a significant fraction of the members become gravi- 
tationally unbound. 

2.3. Numerical Experiments 

In this section we discuss the specific parameters that 
were varied to study the effects of particular initial condi- 
tions on the kinematic structure observed in the clusters. 
A more detailed discussion of each experiment's initial 
set-up and result is reserved until fj3l 

Cluster Geometry and Virial Ratio: We compare the 
evolution of clusters under various assumptions of initial 
geometry and virial balance. Specifically, we compare 
centrally concentrated {p ^ r^^) spherical clusters with 
subvirial {Q — 0.04) and virial {Q ~ 0.5) initial veloci- 
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TABLE 1 
Stellar Distributions 



Description Density Profile Parameters 

Spherical p{^) ={^"'''' ^tl'^ i = r/ro 

< m < 1 m2 = (x/af + {y/bf + {z/cf 
m > 1 a = b < c 

0<m<l m2 = (x/a)2 + (j//b)2 + (z/c)2 
a = b < c 



Uniform Spheroid - 

1/m Spheroid p{m) = <! fj'" ^"^"^ 




TABLE 2 
Embedding Gas Distributions 



Description Density Profile Parameters 

C = r/ro 
{x/a)^ + {y/bf + (z/c)^ 

{x/af + {y/bf + {z/cf 
a = b < c 



Spherical p{0 = Po/^{l + if, 0<^ 

Uniform Spheroid ~ o''' m~>'2 ~ ^ 

1/m Spheroid p(m) = Po/"i, < m 



tics to uniform density elongated clusters with similarly 
subvirial and virial velocities. The spherical clusters arc 
2 pc in radius, and the elongated clusters have axis pa- 
rameters a = 6 = 2 pc and c = 4 pc. Each of these 
clusters had stellar membership N = 2000 and an SFE 
of 50%. We find that only clusters with subvirial initial 
velocities and elongated geometries produce significant 
gradients in the radial velocity along the length of the 
cluster. The observed radial velocity gradients are thus 
a combined effect of global collapse and the projection of 
a non-spherical cluster. Hereafter, the term 'kinematic 
signature' refers to this radial velocity gradient. This 
kinematic signature is discussed in detail in ii3.ll 

With the requirements of subvirial velocities and elon- 
gated stellar/gas distributions identified as prerequisites 
for the kinematic signature, we proceed to complete a se- 
ries of numerical experiments to determine how changing 
cluster parameters in subvirial elongated clusters changes 
the observed structure of the signature. 

Initial Virial Ratio: To characterize the effect of sub- 
virial velocities on the kinematic signature, the evolution 
of centrally concentrated elongated clusters with virial 
ratio Q ranging from 0.04 to 0.15 are compared. These 
virial ratios correspond to average initial velocities that 
are roughly one third to one half of the virial velocity 
and are comparable to pre-stellar clump velocities ob- 
served many s tar forming regions including NGC 2264 
(jPeretto et al.| [2006l and Perseus (jKirk et al.ll2007t ). 

Cluster Elongation: The effect of cluster elongation on 
the evolved velocity structure of the cluster is studied by 
comparing subvirial clusters that range from spherical to 
elongated with aspect ratios ranging from 1 to 4. Recent 
results from the Spitzer Young Cluster Survey indicate 
this range of aspect ratios is appropriate, as the clusters 
in the survey had aspect ratios between 0.53 and 3.88 
(Gutermuth et al. 2008, in preparation). 

Initial Density Distribution: Another experiment com- 
pares subvirial elongated clusters with uniform density 
distributions to those with p ^ to study the effect 
of the density distribution on the strength of and evolu- 



tion of the kinematic signature. 

Star Formation Efficiency: To determine the effect of 
SFE on the kinematic signatures observed in the cluster, 
we compare subvirial elongated clusters with efficiencies 
ranging from 17 to 50%. 

3. SIMULATION RESULTS 

3.1. Radial Velocity Structure Due to Global Collapse 
and Elongation 

We find that radial velocity gradients along the length 
of the clusters are produced by a combination of two 
effects: [1] projection of an elongated or non-spherical 
cluster and [2] subvirial initial velocities. The separation 
of these two effects is nontrivial, as discussed below. 

A cluster that is initially subvirial will collapse as stars 
in the outer parts of the cluster fall toward the clus- 
ter's center of mass. This collapse takes place because 
cluster members that are seeded with subvirial veloci- 
ties are not moving fast enough to remain in orbit (at 
their starting radial positions) around the cluster's cen- 
ter. Instead they fall through the gravitational potential, 
gain kinetic and potential energy, and eventually reach 
an equilibrium state in which the virial theorem is sat- 
isfied. Thus, during the first crossing time, a subvirial 
cluster will collapse significantly as the stellar velocities 
increase. Figure [1] compares the mean cluster radius and 
velocity dispersion as a function of time for spherical clus- 
ters with both virial and subvirial initial conditions. The 
overall cluster collapse and velocity enhancement in sub- 
virial clusters is clearly demonstrated by these plots. 

During the initial collapse phase, the velocity vectors 
of the cluster members are preferentially directed toward 
the center of mass until the stars pass close to the clus- 
ter's center and continue on (mostly) radial trajectories 
outwards. During this first crossing time, some close en- 
counters may cause the individual stellar m otions to devi- 
ate fro m this general pattern. As shown in I Adams et al.l 
(|2006[ ). however, close encounters are relatively rare on 
these short time scales (a few Myr or less) and the aggre- 
gate dynamics are thus dominated by this initial collapse 
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Fig. 1. — The top panel displays the evolution of the mean radius 
of a spherical cluster with virial (dashed curve) and subvirial (solid 
curve) initial velocities. The bottom panel displays the evolution 
of the velocity dispersion of a spherical cluster with virial (dashed 
curve) and subvirial (solid curve) initial velocities. 
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Fig. 2. — Diagram represents an elongated collapsing cluster 
which is viewed along a line of sight not coincident with a pri- 
mary axis. During the initial (collapse-like) phase of a cluster with 
subvirial initial conditions, the upper portion of the cluster will ap- 
pear rcd-shiftcd while the lower portion of the cluster appears more 
blue-shifted. This asymmetry results in a radial velocity gradient 
along the extent of the cluster. 

and re-expansion. 

Elongated subvirial clusters display a gradient in the 
radial velocity along the direction of elongation due to 
this initial collapse. For example, consider a cluster with 
its major axis in the z axis direction and let the system 
be observed along a line of sight that is less than 90° 
from the major axis. For clarity, we define North and 
South to be the projected positive and negative z axes, 
respectively. During the first half of a crossing time, the 
northern part of the cluster appears to be red-shifted 
away from the observer while the southern hemisphere 
is blue-shifted toward the observer. The collapse of an 
elongated subvirial cluster naturally results in a North- 
South gradient in the radial velocities along the length of 
the cluster (see Fig. [2|). The magnitude and direction of 
the radial velocity gradient depends on the line of sight 
chosen to 'observe' the simulated cluster (see §3.3p . 

This kinematic signature in the radial velocities is 
present only if an elongated cluster is viewed along a 
line of sight that is not coincident with a principal axis 
and if the cluster has subvirial velocities. When viewed 
along a minor axis, a subvirial elongated cluster displays 



no North-South gradient in the radial velocities, as the 
stars are preferentially moving perpendicular to the line 
of sight. Viewed along the major axis, the portion of 
the cluster moving away from the observer is projected 
directly onto the region of the cluster moving toward the 
observer, and hence a gradient is not observed. 

Likewise, a virial cluster does not display a strong ra- 
dial velocity gradient, as it is not globally collapsing and 
the stellar velocities have no spatial correlation. A very 
slight gradient is observed in the virial clusters at early 
times; however, it is approximately 10 times weaker than 
that observed in the subvirial clusters and is associated 
with the slight contraction of the simulated virial clus- 
ters as the stellar velocities are redistributed from the 
random isotropic distribution to a true virial distribu- 
tion (see Fig. [1]). 

A cluster that is undergoing global expansion can also 
produce a gradient in the radial velocities similar to that 
observed in the collapsing elongated cluster (Figure [3l 
panel(d)). However, the signature will differ in shape 
because the stars that populate outermost spatial bins 
will not have bound orbits. The stars will not be at the 
furthest point in their orbits and therefore will not have 
reduced velocities as they change direction and return 
to the cluster's center. Therefore, instead of having an 
S-shaped kinematic signature, an expanding cluster will 
have a relatively flat gradient, with no turnover at the 
end points. We also note that rotation of a virial cluster 
may produce a kinematic signature that is similar to the 
one observed in subvirial clusters, but only if the rotation 
takes place around an axis perpendicular to the line of 
sight. 

Figure[3]illustrates the radial velocity signature for var- 
ious projections of spherical and elongated subvirial and 
virial clusters. In each panel the data points correspond 
to the mean radial velocity of the stars in 0.2 pc bins and 
the length of the error bars correspond to the standard 
deviation of the velocities within each bin. For complete- 
ness, we note that this radial velocity gradient calculated 
from the simulated clusters is not sensitive to the choice 
of bin size for bin sizes ranging from 0.05 to 1.0 pc. In 
individual clusters (such as the ONC), the lower limit of 
the bin size is determined by small number statistics. 

It is clear from this figure that the radial velocity signa- 
ture is created by a combination of two effects: projection 
of an elongated cluster and subvirial initial velocities (the 
latter implies global collapse). The bottom right panel 
shows the only cluster that demonstrates this kinematic 
signature. This cluster is elongated, is not viewed along 
a primary axis, and is seeded with subvirial initial veloc- 
ities. 

One way to characterize the strength of a radial veloc- 
ity gradient is by using the amplitude of the cumulative 
velocity distribution, Acrv- The cumulative radial ve- 
locity distribution is created by sorting the individual ra- 
dial velocity measurements by the z' coordinate and then 
producing the cumulative distribution from this sorted 
data set. We normalize this measure by both the number 
of stars in the cluster, N, and the velocity dispersion of 
the cluster, ctrv- This procedure results in a (dimension- 
less) normalized cumulative radial velocity distribution 
that can be more meaningfully compared across clusters 
of different sizes and velocity dispersions. An additional 
advantage of the cumulative velocity distribution is that 
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Fig. 3. — Radial velocity as a function of z' position for a dif- 
fering geometries, initial velocity distributions, and projections for 
simulated embedded clusters. The mean radial velocity (averaged 
over 0.2 pc bins) is in indicated by the data points. The error bars 
correspond to the standard deviation of the radial velocity distri- 
bution within the bins. The panels correspond to differing initial 
geometries, initial velocity distribution, and projections of simu- 
lated embedded clusters: (a) spherical cluster, with subvirial initial 
conditions, (b) elongated cluster, with subvirial initial conditions, 
and viewing angle of 0° from major axis, (c) elongated cluster, 
with virial initial conditions, and viewing angle of 30° from major 
axis, and (d) elongated cluster, with subvirial initial conditions, 
and viewing angle of 30° from major axis. The first distribution is 
shown at time t = 1.5 Myr after the start, whereas the other three 
cases are shown at t = 2.1 Myr. 

it is created from individual stellar radial velocities and 
thus is insensitive to the choice of bin size. 

A cluster with no radial velocity gradient will have (on 
average) the same number of stars with positive and neg- 
ative radial velocities with respect to the cluster's center 
of mass or mean radial velocity. Therefore, the cumula- 
tive velocity distribution in such a cluster will fluctuate 
around zero and be relatively flat as a function of z' . 

On the other hand, a cluster with a strong gradient 
along z' will have preferentially more blue-shifted stars 
at the negative z' end of the cluster center and more 
red-shifted stars at the positive z' end of the cluster. 
Therefore, the cumulative radial velocity distribution will 
be a decreasing function for z' < and an increasing 
function for z' > 0, resulting in a dip in the cumulative 
radial velocity distribution. 

The normalized cumulative radial velocity distribu- 
tions for the clusters depicted in Figure [3] are displayed 
in Figure [H The ranges on the y axes arc held constant 
to emphasize how this distribution varies for subvirial 
and virial initial conditions. Although the virial elon- 
gated cluster in panel (c) of Figure [3] appears to have a 
slight gradient in the radial velocity, the cumulative dis- 
tribution indicates that the strength of this signature is 
more than 7 times weaker than in the subvirial elongated 
cluster. The cumulative distribution for the virial clus- 
ter peaks at Acrv ~ 0.04, whereas that of the subvirial 
cluster peaks at Acrv ~ 0.29. Note that the minimum is 
not observed in the subvirial spherical cluster, panel (a), 
which indicates that subvirial velocities alone are insuf- 
ficient to produce this kinematic signature. As a result, 
in general, both subvirial and non-sphcrical initial con- 
ditions are required to observe the kinematic signature. 

The key feature in the cumulative radial velocity dis- 
tribution signature is neither the growth nor the decay 
of the ftmction, but rather the peak produced by both 
a region of growth and then a region of decay. While 
either the growth, or the decay, of the distribution can 
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Fig. 4. — Cumulative radial velocity distribution (normalized by 
cluster size and velocity dispersion) summed along the i' axis of 
the projected cluster. The panels correspond to differing initial ge- 
ometries, initial velocity distribution, and projections of simulated 
embedded clusters: (a) spherical cluster, with subvirial initial con- 
ditions, (b) elongated cluster, with subvirial initial conditions, and 
viewing angle of 0° from major axis, (c) elongated cluster, with 
virial initial conditions, and viewing angle of 30° from major axis, 
and (d) elongated cluster, with subvirial initial conditions, and 
viewing angle of 30° from major axis. The first distribution is 
shown at time t = 1.5 Myr after the start, whereas the other three 
cases are shown at t = 2.1 Myr. 

be mimicked by constructing the distribution in a veloc- 
ity frame that is significantly different from the average 
velocity of the cluster, both cannot be created in the 
same frame. For instance, if the cumulative radial ve- 
locity distribution for a virial cluster is constructed by 
shifting into a velocity frame that is significantly higher 
than the cluster's average velocity, then almost all stars 
contribute negative values to the distribution, and thus 
the distribution decreases as a function of z' . For the 
purpose of this study, we construct all cumulative radial 
velocity distributions by shifting each stellar velocity by 
the median radial velocity. 

We also note that a flat radial velocity gradient, such 
as that produced by a rotating virial cluster, or an ex- 
panding elongated cluster, can produce a peak in the 
cumulative radial velocity distribution similar to the one 
observed in the collapsing elongated cluster. Therefore, 
both the kinematic signatures presented in Figure [3] and 
the integrated distributions presented in Figure [4] should 
be considered when looking for signs of collapse in an 
observed cluster's radial velocity structure. 

3.2. Velocity Dispersions 

Observations of kinematic signatures in young clusters 
can be compared to these simulated signatures only if 
the observed clusters are large enough that the binned 
radial velocity data do not suffer significantly from small 
number statistics. We investigate trends in the total, ra- 
dial, and tangential velocity dispersions in order to allow 
for comparison to smaller and less studied clusters for 
which only estimates for the total velocity dispersions 
are available. 

In both spherical and elongated subvirial clusters, the 
velocity dispersions increase as a function of time dur- 
ing the initial collapse phase and then decrease as the 
cluster evolves toward virial equilibrium. Clusters that 
experience multiple collapse and re-expansion cycles be- 
fore gas removal also display corresponding cycles in their 
velocity dispersions. The dispersions peak before the 
re-expansion phase and bottom out before the collapse 
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phase, with the peak values becoming smaller with each 
successive cycle. Virial clusters have velocity dispersions 
that increase only slightly during the first 1 Myr, the star 
formation epoch, and then remain relatively flat through 
the rest of the embedded stage. These trends are simi- 
lar to the ones depicted in Figure [1] While an elongated 
cluster's total velocity dispersion does not depend on the 
projection, its radial and tangential velocity dispersions 
do. This dependence on projection angle, 9 is discussed 
in gSXl 

3.3. Kinematic Signatures 

As discussed in t^3.11 only clusters with subvirial ini- 
tial velocities and elongated geometries produce radial 
velocity gradients. The kinematic signature evolves over 
time. During the first 1 Myr, the radial velocity gradient 
arises as stars become free to move through the gravi- 
tational potential of the cluster. The gradient becomes 
larger during the first 'free fall time' while the stars fall 
toward the cluster center. As a majority of the stars 
pass through the cluster center, the radial velocity gra- 
dient decreases and eventually changes sign as the clus- 
ter re-expands. Clusters with sufficiently small crossing 
times will undergo several gradient sign changes before 
gas removal at t = 5 Myr as the cluster size oscillates 
around its equilibrium size. In addition, as the cluster 
evolves toward equilibrium, the strength of the kinematic 
signature decays. 

At t = 5 Myr, the embedding gas potential is removed 
from the cluster which then continues to evolve due to 
gravitational interactions between the stars. Depending 
on the specific kinematic status of the cluster at the time 
of gas removal, the signature considered here may or may 
not remain in the cluster. If the gas is removed while a 
significant portion of the stars are in the re-expanding 
phase, the kinematic signature is amplified and remains 
in the cluster even after 5 Myr of evolution. It is impor- 
tant to note that the details of the gas removal process 
may significantly affect the characteristics of the signa- 
ture, and hence we refrain from making any strong con- 
clusions about the kinematic signature's presence after 
gas removal. 

3.3.1. Effect of Initial Virial Ratio 

We compare centrally concentrated elongated clusters 
(p ^ m~^) clusters with virial parameters Q = 0.04 and 
0.15 to determine how the departure from virial equilib- 
rium effects the strength of the kinematic signature. The 
clusters have N = 2000 members, a SFE of 50%, and axis 
parameters a = 6 = 2 pc and c = 4 pc. We find that as 
the initial virial ratio (and the average starting stellar 
velocity) decreases, the strength of the kinematic signa- 
ture increases. Therefore, the larger the departure from 
virial equilibrium corresponds to a larger kinematic sig- 
nature. We measure the strength of the signature by the 
amplitude of the normalized cumulative radial velocity 
distribution Acrv^ and find that the Q = 0.04 cluster 
had a maximum at Acrv ^ 0.33 whereas the Q = 0.15 
cluster had a maximum at Acrv ^ 0.20. 

In addition, even though the Q — 0.15 cluster had 
larger initial velocities, the Q = 0.04 cluster has a larger 
average velocity dispersion over the embedded stage. 
Specifically, the average velocity dispersions were 0.51 
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Fig. 5. — Amplitude of the normalized cumulative radial velocity 
distribution AcRV, as a function of projection angle d. The data 
arc shown for an elongated cluster with axis ratios (1:1:2), 
uniform stellar and gas densities and subvirial initial velocities at 
the time t = 2.1 Myr. 

km s~i and 0.47 km s"^ for the Q = 0.04 and 0.15 clus- 
ters respectively. 

3.3.2. Effect of Projection Angle 

For a subvirial elongated cluster, the strength of the 
radial velocity signature varies most strongly with pro- 
jection angle. This trend is due to two competing effects. 
First, the most significant collapse occurs along the ma- 
jor axis of the cluster, the z axis. The component of the 
velocity that is observed along the line of sight varies as 
the cosine of the projection angle between the z axis 
and the line of sight. Therefore, the smaller the pro- 
jection angle, the stronger the signature. Secondly, the 
radial velocity signature is only observed when there are 
significantly more red-shifted stars in a declination bin 
than blue-shifted stars (or vice versa). Therefore, a larger 
projection angle will cause the red- and blue-shifted 
populations to be more spatially separated, whereas a 
small projection angle will result in a projected cluster 
whose red-shifted and blue-shifted populations appear to 
overlap. 

Figure [5] displays the amplitude of the normalized cu- 
mulative radial velocity distribution Acrv , as a function 
of projection angle 6, for the uniform density elongated 
cluster with a = 6 = 2 pc and c = 4 pc at t = 2.1 Myr. 
A projection angle between 25 and 35 degrees from the 
z axis produces the strongest radial velocity signature. 
The peak of the distribution is relatively broad, however, 
with a full width at half maximum of approximately 50 
degrees. 

During collapse, the dispersion in the radial velocities 
urv is a decreasing function of the projection angle 9 
from the z axis, whereas the velocity dispersion along 
the z' axis ^w.z' is an increasing function of 9. This result 
occurs because the dispersions are geometrically related 
to each other by (Tflv (6*) = cr||^z/ (90 — 0). These trends are 
representative of the fact that during the collapse phase, 
the velocity dispersions are largest along the principal 
axis. 

In contrast, the plane-of-sky velocity dispersion per- 
pendicular to the z' axis does not depend on the pro- 
jection angle. This result is due to the x, y symmetry 
of the prolate spheroid. Figure [6] displays (jRy (top 
panel), cry j,/ (middle panel), and a^^z' (bottom panel) 
for the same cluster shown in Figure [5l Each curve 
corresponds to a different choice of projection angle, 
6* = 0, 15, 30, 45, 60, 75, and 90 degrees. The relation- 
ship between urv and cry, 2' described above is apparent 
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Fig. 6. — The evolution of the velocity dispersions in an elongated 
subvirial cluster during the embedded stage for various projection 
angles. The top panel displays the radial velocity dispersion ajiy 
as a function of time for 6 = 0, 15, 30, 45, 60, 75, and 90 degrees. 
The middle panel displays the velocity dispersion parallel to the 
z' axis, ,/ for the same projection angles. Note that a-RYiS) = 
cry 2/(90 — 0) by definition. The bottom panel shows the velocity 
dispersion in the plane-of-sky perpendicular to the z' axis, a-± ;.i. 
On account of the x, y symmetry of a prolate spheroid, crj_^-/ does 
not vary as a function of 8. The data presented here are from a 
uniform spheroid cluster with subvirial initial velocities and axis 
parameters a = fe = 2 pc and c = 4 pc. 

in Figure [6l 

3.3.3. Effect of Cluster Elongation 

To observe differences caused by the amount of elon- 
gation, we compare two subvirial clusters with differing 
axis ratios. Specifically, wc compare a cluster with an as- 
pect ratio of 0.5 (a = 6 = 2 pc, c = 4 pc) to one with an 
aspect ratio of 0.25 (a = 6 = 1 pc, c = 4 pc). Each of the 
clusters has 2700 members, a density profile p ~ m~^, a 
star formation efficiency of 33%, and a virial parameter 
Q = 0.04. 

We find that the more elongated cluster evolves on 
a slightly shorter time scale and that the mean radial 
velocities in a given declination bin are larger. As a re- 
sult, the kinematic signature discussed in ^S.ll is stronger. 
For comparison, the maximum mean radial velocity in a 
declination bin is approximately 1 km s~^ for the less 
elongated cluster and 2 km s~^ for the more strongly 
elongated cluster. Therefore the overall change in radial 
velocity is approximately 4 km s~^ over the length of the 
cluster in the strongly elongated geometry. The average 
velocity dispersion is also larger in the more elongated 
cluster, a = 1.60 km s~^, compared to the less elongated 
cluster where a = 1.18 km s~^. 

3.3.4. Effect of Initial Density Distribution 

The kinematic signatures also vary as a function of 
initial density distribution. We investigate this varia- 
tion by comparing subvirial {Q = 0.04) elongated cluster 
with uniform density distributions to those with centrally 
concentrated p ~ density profiles. Each cluster con- 
tains N — 2000 members, has an SFE of 50%, and has 
axis parameters a = & = 2 pc and c = 4 pc. 

Elongated clusters that are more centrally concen- 
trated have central regions which evolve on shorter time 



scales than the outer regions. This trend can be seen by 
considering the gravitational contraction of an centrally 
concentrated sphere with radius tq and density profile 
p ~ . The collapse is inside-out because free fall 
time for a particle initially at a position r is given by 
tff ^ ~ ^/r, and so the centr al regions of the sys- 

tem e volve on shorter timescales. iGoodman &: BinnevI 
(|1983[ ) showed that this general argument holds true for 
spheroidal clusters as well. 

On account of this inside-out collapse, the center of 
a subvirial cluster will collapse and begin to re-expand 
while the outskirts of the cluster are still collapsing. As 
a result, the red and blue shifted populations will not be 
as spatially separated as in the uniform density cluster 
shown in Figure [31 Panel (d). Instead, the generally red- 
shifted z' > portion of the cluster will contain many 
blue-shifted stars near the center, and this addition acts 
to blue-shift the average radial velocity in the central 
z' bins. On the other side of the cluster's center, red- 
shifting may occur for the same reason, resulting in a 
kinematic signature that is flat or even oppositely sloped 
(with respect to the general trend) near z' = 0. 

In addition, the inside-out collapse causes the aver- 
age velocity dispersion during the embedded stage to be 
higher in centrally concentrated clusters than in uniform 
density clusters. This effect occurs because the peak in 
the velocity dispersion, which occurs during the cluster's 
collapse and re-expansion, is broader (in time) in the 
p ~ clusters than in the p ~ po clusters. In a cen- 
trally concentrated spheroid, the stellar members reach 
their highest speeds at a time that depends on their ini- 
tial position within the cluster. In contrast, in a uni- 
form spheroid, the stellar members reach their highest 
speeds at roughly the same time. Therefore, the veloc- 
ity dispersion in a centrally concentrated cluster remains 
higher for a longer period of time. This finding has ob- 
servational consequences: initially centrally concentrated 
clusters are more likely to be seen with high velocity dis- 
persions. 

3.3.5. Effect of Star Formation Efficiency 

To investigate the effect of SFE on the kinematic signa- 
tures observed, we compare elongated clusters with SFE 
of 17%, 33%, and 50%. Each cluster has a p ~ to~^ stel- 
lar density profile, axis parameters a = 6 = 2 pc and 
c — 4 pc, and subvirial Q ~ 0.04 initial velocities. The 
clusters with SFEs of 17% and 13% have N = 2700 stars, 
whereas the cluster with a 50% SFE has N ~ 2000 stars. 

Clusters with lower star formation efficiencies evolve on 
shorter timescales due to the higher velocities associated 
with the deeper potential wells of the gas component. In 
addition, both the amplitude Ackv and the total veloc- 
ity dispersion increase as SFE decreases. 

3.3.6. Effect of Extinction due to High Column Density 

Most very young clusters are still embedded in the 
molecular cloud from which they formed. In cases where 
optical spectra are used to determine stellar kinemat- 
ics, extinction from the embedding clouds can poten- 
tially change the observed kinematic signatures in these 
clusters. For example, consider a cluster that is roughly 
shaped like a prolate spheroid with its North end tipped 
toward the observer. The Southern region is farther 
away, and is collapsing toward the center of mass, and 
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Fig. 7. — Same as Figure [2l but indicates portion of cluster that 
is assumed to be obscured due to extinction. 

thus appears blue-shifted with respect to the stars in the 
northern region. If even a modest fraction of the cluster 
members that are farthest from the observer are not in- 
cluded in the observed data set due to extinction by the 
embedding cloud, this blue-shifted population would ap- 
pear less blue-shifted than if it were not obscured. These 
observed selection effects may be mitigated by using H 
or K band spectra to select targets. 

As an example, consider an embedding cloud that has 
a dust opacity of Ky = 200 cm g^^ at optical wavelengths 
and a mean molecular mass fj, = 2.4m//. In order to have 
an optical depth with the value of unity, the minimum 
column density is required to be Ncoi ~ 1-28 x 10^^ cm^^. 
A typical molecular cloud has number densities ranging 
from n = 10^ — 10^ cm^'^. Therefore, if the sources are 
being observed through more than 0.4 — 4 pc of molecular 
cloud, they would be undetected and hence removed from 
the sample. 

To test this hypothesis, we re-analyzed the simulated 
cluster data by omitting stars that are behind an "ob- 
scuring plane" (see Fig. [7|) which includes a column den- 
sity threshold. In this model, the obscuring plane is in- 
troduced with an orientation normal to the line of sight 
and at a distance < dobs < 5 pc beyond the center 
of the cluster from the observer. We found that if only 
a modest number of stars in the Southern region of the 
cluster are removed from the sample (due to the fact 
that they are beyond the obscuration plane) and their 
velocities are not included in the velocity versus posi- 
tion plot, the Southern portion of the cluster will have a 
significantly less steep velocity gradient compared to the 
Northern region. Specifically, if 10 — 15% of the stellar 
population is unobservable, the kinematic signature in 
the southern half of the simulated cluster is completely 
washed out. Figure [8] displays the radial velocity sig- 
natures for a cluster with various amounts of the stars 
extincted. In the top panel, dobs = 0.6 pc and 32% of 
the cluster members are removed from the sample. The 
radial velocity signature in the middle panel is observed 
when dobs = 1-4 pc and 12% of the cluster members are 
thus obscured. For comparison, the bottom panel is the 
radial velocity signature produced when all stars are in- 
cluded in the sample. 

4. COMPARISON TO OBSERVATIONS 
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Fig. 8. — Each panel displays the radial velocity as a function 
of z' position in a uniform density subvirial elongated cluster for 
different locations of the obscuring plane. The mean radial veloc- 
ity (averaged over 0.2 pc bins) is in indicated by the data points. 
The error bars correspond to the standard deviation of the radial 
velocity distribution within the bins. In the top panel, dobs = 0.6 
pc and 32% of the cluster members are removed from the sample. 
The radial velocity signature in the middle panel is observed when 
fi-oba = 1-4 pc and 12% of the cluster members are thus obscured. 
The bottom panel is the radial velocity signature produced when 
all stars are included in the sample. Even modest amounts of ex- 
tinction can wash out the structure of the kinematic signature. 

Recent kinematic studies of the Orion star-forming re- 
gion have determined radial velocities for a large sample 
(~ 1200) of visible sources within the region. In this sec- 
tion we discuss the kinematic results of these studies in 
light of the numerical simulations presented above. 

The ONC region is a good environment in which to 
look for kinematic signatures such as the ones observed 
in our simulated clusters for many reasons. First, the 
ONC has a large population and is close enough to be 
well studied so that many of the cluster members have 
measured radial velocities. The large sample size of the 
cluster allows the data to be binned in declination while 
maintaining a reasonable number of data points in each 
bin, so that the results suffer only mildly from small num- 
ber statistics. Furthermore, the ONC is visibly elongated 
in projection (jHillenbrand fc Hartmannlll998f ). which is 
one of the requirements for the kinematic signature to be 
observed. 

Finally, although opinions vary, there are significant 
arguments supporting the assertion that the ONC re- 
gion as a whole is young and is estimated to be less 
th an one crossing tirn e old. Observations of the ONC 
bv iFiiresz et al.l ()2008D and more recently by T09 iden- 
tify spatially coherent kinematic structure in the stellar 
distribution that closely matches the observed kinematic 
str ucture of the gas in the region as measured in ^^CO 
bv lBallv et al.l (jl987f ). The authors argue that this cor- 
relation indicates the cluster is not dynamically relaxed 
and that in fact, the region is less than a crossing time 
old. The results of the previous section indicate that 
the strength of the kinematic signature peaks before the 
cluster is a crossing time old. Therefore, the strongest 
observed signature will occur in clusters that are less than 
a crossing time old, such as the ONC. 
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with radial velocities measured by IIFuresz et al.ll2"008l . T09) that 
are spatially associated with the filament of dense gas. In the 
right panel, the median radial velocity of the stellar population 
is plotted for each declination bin. Bin width is indicated by the 
vertical extent of the bars. The horizontal error bars indicate the 
velocity dispersion in each bin, where the dispersion is defined by 
the distribution's full width at half maximum divided by 2V21n2. 
There is a noticeable radial velocity gradient in the Northern ONC, 
which is qualitatively similar to the one observed in the simulated 
subvirial clusters. 

4.1. Observations 

The radial velocities cited in this study were deter- 
mined from multiple observations using Hcctochellc on 
MMT and MIKE fibers at Magellan The Hectochelle 
observations from iFiiresz et al.l (|2008f l have been com- 
bined with additional Hectochelle and MIKE observa- 
tions (T09) to produce a fist of ~ 1200 sources with 
radial velocity measurements in the region surrounding 
the ONC. A detailed discussion of the observations, data 
reduction, spectral fitting, and radial velocity measure- 
ments is provided in the observational papers. In addi- 
tion, much care has been taken to identify binaries and 
non-cluster members from the ONC data, and this pro- 
cedure is detailed in T09. 

4.2. Kinematic Signatures 

As discussed in these observational papers 
(|Furesz et al.l l2008l T09), the stars and gas in the 
region surrounding the ONC show similar North-South 
velocity gradients. We consider the ONC members that 
are near the molecular cloud filament by selecting the 
RA range 84.0 to 83.5, and remove identified binaries 
from the sample. The distribution of sources is shown in 
the left panel of Figure [9l In the right panel of Figure 
[S] the binned median velocities are plotted as a function 
declination. To calculate the median velocity in each 
declination bin, a histogram of the radial velocities is 
created. The median is calculated ignoring all bins 
that have less than half of the maximum value in the 
distribution. The medians calculated without localizing 
on the filament differ from those shown here by a few 
0.5 km/s shifts to the red or blue. 

Figure [9] is analogous to Figure [3] for the simulated 
cluster data. The North-South velocity gradient is clear, 
though the gradient is steeper in the northern part of 
the ONC (north of declination ^ —5.5) than it is in the 
southern region. The observed kinematic structure in 
the ONC region is qualitatively similar to the kinematic 
structure of elongated subvirial clusters viewed off axis 
(bottom right panel of Fig. [S] 
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Fig. 10. — Normalized cumulative radial velocity distribution 
summed along declination in the Orion star forming region. This 
cumulative distribution is created from the ONC members with 
radial velocities within 3 cr of the m edian radial velocity of the dis- 
tribution. From IFurcsz et al.l 1)20081 ) , the radial velocity dispersion 
in this region of the ONC is cr = 3.1 km s~^. The cumulative 
radial velocity distributions for the stars north of declination —6 
and —5.5 are indicated by the dashed lines. 

Figure [10] presents the cumulative radial velocity dis- 
tribution observed in the ONC (analogous to Figure [3]) 
for the entire sample (solid curve) as well as the distri- 
bution for only those stars north of declination —6 and 
—5.5 (dashed curves). This distribution is composed of 
the radial velocities of the ONC members that have ve- 
locities within 3 fj of the median radial velocity of the 
dist ribution of the me dian velocity, where cr = 3.1 km 
s~^ (jFiiresz et al.ll200"8l ). The constrained data set (with- 
out the declination cuts) includes approximately 89% of 
the sources shown in the left panel of Figure [9l To pro- 
duce these plots, the radial velocities were shifted by the 
median observed velocity and normalized by the velocity 
dispersion of and number of stars in the data set. This 
shift and normalization allows the observed data to be 
directly compared to the cumulative radial velocity dis- 
tribution in the simulated clusters, where the velocities 
are measured in the center of mass reference frame. The 
kinematic signature peak observed in our simulated clus- 
ters is also observed in the ONC, over the extent of the 
cluster, as well as over just the northern portion. The 
ONC data is significantly less smooth than the theoreti- 
cal results. This difference is mostly due to the summing 
of 50 simulations used to produce the theoretical results 
shown in Figure [H 

4.3. Comparison to Simulation Results 

The kinematic structure observed in the ONC may be 
understood in light of the simulated cluster kinematics. 
The correlation between the motion of the stars and the 
gas in the ONC indicates that the region is fairly young 
and that the stars and gas are still collapsing. The ob- 
served radial velocity gradient also supports the theory 
that the ONC is dynamically young and that the region 
is less than a crossing time old. Also, the signature sug- 
gests the stars in the region are formed with subvirial ve- 
locities, an initial condition required by the simulations 
to produce the structure observed. The ONC region is 
elongated and the combination of elongation, subvirial 
initial velocities, and extremely young dynamical age can 
account for the kinematic structure observed as demon- 
strated by the simulations presented in this study. This 
interpretation requires that the Northern region of the 
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ONC be closer to the observer than the Southern region. 
However, at this time, the dista nce to the O NC as a 
whole is still poorly constrained ([Jeffries! [20071 . and ref- 
erences therein) and determining distance differences on 
the order of 10 pc between the Northern and Southern 
regions of the ONC is not yet possible. 

The kinematic signatures observed in the ONC are 
qualitatively similar to those observed in the elongated 
subvirial clusters. The simulated cluster that is phys- 
ically most similar to the ONC region is the centrally 
concentrated elongated cluster with axis ratio: 2:2:8 
pc and a low star-formation efficiency of ^ 17%. This 
cluster displays a kinematic signature similar in magni- 
tude to the observed signature in the ONC. Specifically, 
at t = 0.7 Myr, the simulated cluster has a radial ve- 
locity gradient of ^ 2 km s~^ pc~^, and Acrv ~ 0.15. 
In addition, the simulated cluster has a radial velocity 
dispersion of ~ 3 km s~^. In comparison, the Orion 
data have a radial velocity gradient of ~ 5 km s^^ over 
0.6° in declination; at a distance of 420 pc, the param- 
eters are approximately a gradient of 1.1 km pc~^, 
-^CBV ^ 0.2 and a radial velocity dispersion of 3.1 km 
s-i.' 

Proper motion studies in Orio n (such as the one con- 
ducted bv lJones fc Walkerlll988[ ) found plane-of-sky ve- 
locity dispersions along the cluster's major axis that are 
somewhat larger than the dispersions perpendicular to 
the major axis. In the notation adopted in this paper, 
= 2.63 ± 0.9 km s"^ and aj_^^, = 2.03 ± 0.11 km 
s~^. iFiiresz et al.l ()2008l ) determined a radial velocity 
dispersion of auv = 3.1 km s~^. This relatively large 
difference between any and (7±_z' is indicative of some 
type of global collapse because virial elongated clusters 
have at most modest differences between any and <7±^z'- 
In our simulated virial clusters 1 < <7rv /o'±,z' < 1-3, 
whereas in the subvirial clusters 1 < (Jrv /o'±^z' < 2.3. 

The ONC region displays significant mass segregation 
among the most massive stars (to > 2Mq) which reside 
in the Trapezium at the center of the ONC. Previous 
numerical studies of the ONC cluster indicate that the 
mass segregation must be primordial as the cluster is not 
old enough for dynamical mass segregation to account fo r 
the presence of the Trapezium ([Bonnell fc Davieill998[ ). 
Our simulated clusters begin with somewhat different ini- 
tial conditions, by focusing on subvirial initial velocities 
and cluster elongation, but arrive at the same conclusion: 
the region is not old enough for dynamical mass segre- 
gation to have taken place. In the simulated clusters, we 
compared the radii that enclosed 25, 50, and 75% of the 
stellar mass in different mass bins. The results of these 
simulations showed no difference in the radii of different 
mass bins during the first 1 — 3 Myr of evolution, indi- 
cating that dynamical mass segregation had not yet oc- 
curred. Subvirial initial velocities will result in a higher 
central density at earlier times on account of global col- 
lapse, but even with an increased density and thus higher 
interaction rate, the ONC is still too young for significant 
mass segregation to have occurred. 

In addition to the lack of mass segregation in the spa- 
tial distribution of stars within the cluster, no significant 
differences in the kinematic distributions are observed 
as a function of stellar mass in the simulations. The 
masses of stars in the Orion region are not well enough 



constrained to compare this result of the simulations to 
the observations at this time. However, future studies 
of young stars in this region may yield mass information 
that, combined with kinematic data, will be able to test 
this prediction. 

As seen in Figure O the radial velocity gradient ap- 
parent in the northern region is not as evident in the 
region south of the ONC. The asymmetry may indicate 
that the center of the gravitational collapse is not the 
Trapezium, but rather is located slightly north of the 
Trapezium. The collection of cumulative radial velocity 
distributions displayed in Figure [TOl however, indicate 
that the signature is present even if the center of the 
gravitational collapse is assumed to be somewhat north 
of the Trapezium. The asymmetry could also be due to 
a north-south asymmetry in the initial conditions that is 
not well represented by the idealized simulations of fJ31 
The star-forming region north of the ONC is located at 
the edge of the integral-shaped molecular cloud, whereas 
the region south of the ONC extends into a larger com- 
plex of molecular clouds. Therefore the dynamics in the 
northern region is better suited to our approximate of 
an isolated elongated centrally concentrated cluster than 
those of the south. 

5. CONCLUSION 

This study presents A^-body simulations that explore 
the kinematic signatures produced by asphericity and 
subvirial initial velocities in young embedded clusters. 
We have identified a robust kinematic signature, in par- 
ticular a gradient in the radial velocities, which is natu- 
rally produced by elongated subvirial clusters. We char- 
acterize the properties of this kinematic signature as a 
function of initial conditions. Specifically, we compare 
the kinematics of evolving spherical and aspherical clus- 
ters and with both virial and subvirial initial conditions. 
We consider changes in the kinematic signature due to 
differing amounts of cluster elongation, star-formation 
efficiency and departure from virial equilibrium, and dis- 
cuss possible sample bias due to extinction in observed 
clusters. Finally we compare the signature displayed in 
our simulated clusters to kinematic data from the ONC. 
The main results of our work are summarized as follows: 

• Elongated clusters with subvirial initial velocities 
display a gradient in the radial velocities as a func- 
tion of projected position along the cluster. Both 
aspherical initial geometries and subvirial initial ve- 
locities are required to produce this kinematic sig- 
nature (see Fig. [3|). 

• The strength of the kinematic signature increases 
during the first free-fall time as the cluster col- 
lapses, and then decreases as the stars pass through 
the center of the cluster. The gradient changes 
sign as the cluster re-expands, and the amplitude 
of the cumulative radial velocity distribution de- 
creases (see Fig. |4|) as the cluster evolves (see also 
Fig. [J). 

• The strength of the kinematic signature varies most 
sensitively as a function of the projection angle, as 
measured from the major axis of the cluster. The 
signature is weakest at projection angles of 0° or 
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90° and strongest at projection angles of ' 
(see Fig. [S|). 



• Making the cluster more elongated increases the 
gradient of the radial velocity across the cluster 
while leaving the velocity dispersions roughly sim- 
ilar. 

• The initial stellar (and gas) density distribution af- 
fects the rate at which the cluster evolves and thus 
the time scale on which the kinematic signature 
evolves. In addition, centrally concentrated clus- 
ters have higher velocity dispersions than uniform 
density clusters. 

• The kinematic signature is sensitive observational 
selection effects. Extinction in an embedded clus- 
ter can preferentially deselect stars further from the 
observer and will thus affect the kinematic signa- 
ture. A modest amount of extinction that removes 
10% of the cluster members will wash out the kine- 
matic signature in the region of extinction (see Fig. 

[Zl). 

• The asymmetric kinematic signature is qualita- 
tively similar to the observed kinematics of the 
Orion star-forming region, suggesting that in ad- 
dition to being elongated (as observed) a signifi- 
cant fraction of the cluster members started with 
subvirial initial velocities (see The large gra- 
dient observed (~ 1 km s~^ pc~^) indicates that 
the cluster is dynamically young. This result is 
consistent with previou s independent claims (e.g., 
iBonneU fc Davieslfl998f l. 

Previous studies have shown that the star-formation 
efficiency of clusters and the initial velocity distribu- 
tion have signific ant effects on the long term evolu- 
tion of a cluster jAdamd l2000t iBoilv fc Kroupal l2003l : 
I Adams et al.ll2006[ ). i.e., the cluster bound fraction and 
stellar interaction rates. This work emphasizes the fact 
that these conditions affect the short term evolution of 
the cluster as well. In addition, this work is signifi- 
cantly different from previous studies in that it considers 
non-spherical initial conditions in both the stellar and 
gaseous components in young embedded clusters. Al- 
though many examples of clusters with non-spherical ge- 
ometries exist and a considerable amount of work has 
considered orbits of individual stars within axisymmetric 
and triaxial potentials, little work has been done theoret- 
ically to understand how these clusters, taken as a whole, 
differ from spherical clusters. We have shown that elon- 
gated clusters naturally produce observable kinematic 
signatures that depend on the initial starting velocities 
of the cluster members and the projection of the cluster 



25 — 35° onto the plane of the sky. These signatures may allow us 
to identify, or at least constrain, the initial conditions in 
star-forming clusters that still retain most of their cm- 
bedding gas. 

The kinematic signatures found in the elongated sub- 
virial cluster simulations may shed some light on kine- 
matics within nearby young clusters. As one example, 
this work shows that the gradient produced in simulated 
clusters is qualitatively similar to that observed in the 
radial velocities of the ONC members. Although the 
ONC kinematic data display some additional asymmetry 
not observed in the simulated clusters, this asymmetry 
is likely to be due to the more complicated environment 
and feedback processes which are not well represented 
by the simulation. Nevertheless, the general structure 
and magnitude of the observed kinematic signature in 
the ONC may be explained theoretically, provided that 
the stars are formed with subvirial initial velocities. 

Current instruments such as the Hectochelle and 
MIKE have helped make large spectroscopic studies of 
many objects more efficient and have produced many ra- 
dial velocity studies of stellar clusters as well as larger 
globular clusters and dwarf galaxies. These instruments 
have provided insight into the kinematics of astronomi- 
cal objects ranging from the rotation of individual stars 
to the large scale dynamics of galaxies. In the next 
decade, the GAIA mission will provide the opportu- 
nity to study the three dimensional kinematics of clus- 
ters with exquisite detail. For example, proper motions 
and radial velocities will be measured for stars in ONC 
brighter than 20th magnitude. This study will result 
in accurate transverse velocities to combine with the ra- 
dial velocities supplemented by ground-based observa- 
tions. As more detailed kinematic data become avail- 
able for the ONC and other young clusters, they should 
be used in conjunction with cluster simulations such as 
those presented here to understand the initial conditions 
necessary to produce the observed kinematics. A more 
detailed understanding of the initial conditions can then 
inform cluster formation and evolution models, with the 
overarching goal of better understanding star and planet 
formation within young stellar clusters. 
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